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■ Before a quantum-mechanical calculation involving electromagnetic interactions 

is performed, a choice must be made of the gauge to be used in expressing the 

>. 

potentials. If the calculation is done exactly, the observable results it predicts will 
be independent of the choice of gauge. However, in most practical calculations 
approximations are made, which can destroy the gauge invariance of the predictions. 
We compare here the results of coupled-channel time-dependent relativistic Coulomb 
excitation calculations, as performed in either Lorentz or Coulomb gauges. We find 
significant differences when the bombarding energy per nucleon is > 2 GeV, which 
indicates that the common practice of relying completely on the Lorentz gauge can 
be dangerous. 



I. INTRODUCTION 

Coulomb excitation is a collison process in which the predominant projectile-target inter- 
action is electromagnetic. The electromagnetic interaction is expressed in the Hamiltonian, 
and therefore also in the Schrodinger equation, in terms of the electromagnetic potentials, 
(<p(v, t), A(r, t)). The potentials are subject to gauge ambiguity, since they are not uniquely 
determined by the charge and current distributions that create the electromagnetic field. 
Equivalently, if a gauge transformation generated by A(r, t) is performed, the new poten- 



2 



tials (<f'(r, t), A'(r, t)) defined by 

^(r.f) = y?M)-i^M ( i.i a) 

A'(r,t) = A(r,t) + VA(r,t) (1.1b) 

yield the same (E(r, t), B(r, t)) as did (y?(r, t), A(r, i)) Q. Thus (y>(r,*), A(r, *)) and 
(y/(r, t), A'(r, t)) are both consistent with Maxwell's equations and the given charge and 
current distributions. Since the Schrodinger equation depends explicitly on the potentials, 
it is changed by a gauge transformation, and so is the wave function which is the solution 
of the Schrodinger equation. However, when an observable quantity (such as a transition 
probability) is calculated from this wave function, the result is invariant under a gauge 
transformation, even though the wave function is not. 

The gauge invariance of an observable applies only if the exact solution of the Schrodinger 
equation is used in its calculation. If an approximate solution is used, then the calculated 
observable quantity may or may not be gauge invariant. A non-gauge invariant prediction 
is a serious drawback for a theory, since there is generally no a priori reason to choose one 
gauge rather than another. 

An example of an approximate solution that is gauge invariant is given by first-order time- 
dependent perturbation theory, 2]. Suppose we want to estimate transition probabilities in 
the target under the influence of a time- dependent external potential V(r,t) provided by 
the projectile. If we use first-order time-dependent perturbation theory, then the predicted 
transition probability from target state <fr a at t = — oo to target state <p 1 at t = oo is given 
by 

f°° dt 

/ -r * V ya (t) 

J-oo IK 

where 

u-y a = (e 7 — e a )/h (1.2a) 
V^t) = <0 7 I V(r,t) \ 4> a > . (1.2b) 

Thus the first-order transition probability is proportional to the square of the modulus of the 
(onshell) Fourier component of V ia (t). If p ia and J 7a (r, t) are the target transition 
charge and current densities associated with the 4> a — > </> 7 transition, the electromagnetic 
interaction V ia {t) is 



p 7Q (r)<y?(r,t) - ^J 7a (r) • A(r,t) 



;i.3) 
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Target charge conservation can be expressed by the relation 



V-J 7a (r,t) + iuj ia p ia (r) = 0. 

This can be used to show that the uj Fourier component of V ya (t), 

dt 



ili 
dt 

ih 



iu>t 



Vy a (t) 



d 3 r 



is changed by the gauge transformation 

(uj - uj ia ) 

he 



dt 



p ia (r)ip(r, t) - ~J 7a (r, t) ■ A(r, t) 
generated by A(r, i) by 
d 3 re iujt p, /a (r)A{r,t). 



(1.5a) 
;i.5b) 



Thus the onshell, to = uj ja , Fourier component V 7a (uj 7a ), which determines the first-order 
transition probability, is unchanged by a gauge transformation. 

There are other approximations which depend only on the onshell Fourier transform 
of the interaction. For example, in the model in which the target is represented by a 
harmonic oscillator and the interaction is assumed to be linear in the oscillating variable 
or its conjugate momentum, an exact solution of the time-dependent Schrodinger equation 



is available 



]. It depends only on the onshell Fourier transform of the time-dependence 



of the interaction. Also, the FWW "method of virtual quanta" (see, e.g., Reference (J]) is 
essentially first-order perturbation theory, in which the perturbing electromagnetic field is 
approximated by a plane wave. In both theories, the gauge independence of the onshell 
Fourier transform of the interaction guarantees that their predictions are gauge invariant. 

However this gauge invariance is not guaranteed if the full electromagnetic interaction 
is used, and if one wishes to go beyond first-order perturbation theory. For example, a 
coupled-channel approach to the time-dependent Schrodinger equation has been used in an 
attempt to get a better description of Coulomb excitation to multi-phonon states of giant 

UJj . One begins by choosing a set of target states which can be 
expected to play a significant role in the reaction, and then solving the Schrodinger equation 
within the Hilbert space defined by these states. This coupled-channel approach depends 
upon the entire cu-dependence of Vy a (uj), not only on its value for uj = lu joi . If all the target 
states were included in the calculation, then the coupled-channel solution of the Schrodinger 
equation would be exact, and calculated observables would be gauge invariant. However if, 
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TABLE I: Comparison of the corresponding properties of the potentials in the Lorentz and Coulomb 
gauges, due to the projectile charge density, p p (r,t) and current density, J p (r, i). 

Lorentz gauge Coulomb gauge 

V-A L + if ^ = V-A c = 

VV " ^|W L = ~^p P VV = -4vr/5 p 



as is always the case, the set of target states is truncated to make the calculation feasible, 
then the result of the coupled-channel calculation is approximate, and therefore calculated 
observables are not gauge invariant. The main concern of this paper is with the lack of gauge 
invariance of the predictions of coupled- channel time-dependent solutions of the Schrodinger 
equation governing relativistic Coulomb excitation. Some related discussions of this subject 
have been given by Baltz, Rhoades- Brown, and Weneser |l2j, Rumrich and Greiner 43] and 
Kobe and Kennedy [L4J|. 

Q 

The classic paper on relativistic Coulomb excitation 2] by A. Winther and K. Alder 
(WA) described the electromagnetic influence of the projectile on the target using classical 
electromagnetic fields in the Lorentz gauge. Since the main concern of WA was with first- 
order time-dependent perturbation theory, their calculated transition probabilities and cross- 
sections were gauge invariant. Subsequent workers in this field who used coupled-channel 
methods continued to use classical fields and the Lorentz gauge. However, as explained 
above, the extent to which the results of these calculations are gauge invariant is not certain. 

When studies are done in which the electromagnetic field is quantized, it is common to 
use the Coulomb (or transverse) gauge. The reason is that the field quanta (photons) will 
then be purely transverse. Table I compares the properties of the potentials in the Lorentz 
and Coulomb gauges. 

We will consider two models for the radial density of the projectile, as seen in its own 
rest frame: a finite spherically-symmetric charge distribution, and a point charge. The 
expressions for the finite charge distribution will be presented as Fourier transforms, which 
are convenient for the multipole expansions needed when the initial and final nuclear states 
are angular momentum eigenstates. The point charge expressions are presented as functions 
of t, which can be easily interpreted. Numerical calculations show that there is very little 
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difference between the predictions of the two models. 

In Section II we will find the generator A(r, t) which will take us from the commonly used 
Lorentz gauge to the Coulomb gauge in the classical-field approach to relativistic Coulomb 
excitation. This is already known for a point projectile [12j, but we will derive the expression 
appropriate to a projectile of finite size. Section III presents the interaction potentials 
calculated in these two gauges. In Section IV we compare multipole expansions of these 
potentials, and show that the Coulomb gauge potential is free of a divergence that appears 
in the Lorentz gauge potential at high bombarding energy. Sections V and VI apply these 
formulae to the excitation of multiphonon states of the giant dipole resonance in 40 Ca, as 
a result of bombardment by 208 Pb nuclei. Section VII presents our conclusions and some 
general observations about gauge invariance. 

II. THE GAUGE TRANSFORMATION CONNECTING THE LORENTZ AND 

COULOMB GAUGES 

We follow the standard approach to relativisitic Coulomb excitation as proposed by WA. 
The projectile nucleus is assumed to travel along a straight-line orbit parallel to the z axis, 
with impact parameter b, at constant speed v. The magnitude of the impact parameter is 
large enough so that nuclear interactions between the target and projectile are negligible. 
Because of the assumed large projectile momentum, the electromagnetic impulse the pro- 
jectile receives due to its interaction with the target has little effect on its trajectory, so the 
projectile maintains its constant speed and impact parameter throughout the collision. As 
the projectile passes, the target nucleus feels the time-dependent projectile electromagnetic 
fields, which induce transitions between the quantum states of the target. 

We seek the gauge function A(r, £) which generates the gauge transformation 



between the potentials satisfying the Lorentz and Coulomb conditions listed in Table 1. 




(2.1b) 



(2.1a) 
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A. Point projectile of charge Zpe. 



For a point projectile of charge Z p e, we have a charge density given by 

pp(r, t) = Z P e 5(r - b - zvt) = Z P e 5 (p - b + z(z - vt)) . (2.2) 

If this charge density is used on the right-hand sides of the equations for tp in the second 
row of Table I, we get the the solutions 



<^ L M) 

<^(r,t) 



Z p e 



+ (*-«*)'■ 



Z p e 



Baltz, Rhoades-Brown and Weneser 

-t 



(2.3a) 



(2.3b) 



121 ] used these potentials and Equation (2.1b) to obtain 
AM) = c / [^(r,t')-^(r,t')]dt' 

J — oo 

1 1 



Z p ec 



(/Q-b) 2 



+ (z — vt 



A 2 



v/(p-b) 2 + (^-i;f) 2 



Z p e- log 



(^-z) + ; /^^ + (z-^)2 



(2.4) 



(vt - z) + vFf+MF 

Using this A(r, t) and Equation (lb), we can calculate the vector potential in Coulomb 
gauge. In contrast to the Lorentz gauge vector potential, it has a component perpendicular 
to the z direction: 

1 1 



(A c 

(a c ; 



Z p e 



Z p e-(p-b) 



y/(p - b) 2 + (vt - z) 2 7v /(p-b) 2 + 7 2 (t;t-z) 2 

1 



(2.5a) 



y/(p-by + ^(vt~zy (<y(vt -z) + v /(p-b) 2 + 7 2 (^-^) 2 ) 



y/( P ~ b) 2 + (vt ~ Zf ((Vt -Z) + y/( P - b) 2 + (vt ~ Z)A 



(2.5b) 



while 



(A L ) 2 
(A L ) ± = 0. 
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v v 

-ip L (x, y, z, t) = Z p e . 

c 1 ; p c v /(p-b) 2 + 7 2 ( U t 



(2.5c) 
(2.5d) 
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It can be verified by direct calculation that this A c (r, t) is solenoidal, V • A c (r, t) = 0, as 

specified in Table I for the Coulomb gauge. 

The presence of 7 in the Lorentz gauge scalar and vector potentials (2.3a, 2.5c) has 
the effect of decreasing the effective duration of the time-dependent impulse experienced 
by the target. Note that the Coulomb gauge potentials (2.3b, 2.5a,b) have terms that are 
independent of 7. As a result, the Coulomb gauge pulse occurs over a longer time interval 
than does the Lorentz gauge pulse, and so the Coulomb gauge pulse is more adiabatic. This 
will be illustrated in Figures 1 and 2 below. 

B. Finite spherically-symmetric projectile of charge Zpe. 

The method used in Section II. A does not work for a finite-sized projectile. Whereas 
y? L (r, t) is still given by the Equation (2.3a) (outside the projectile), the expression for 
y? c (r, t) is more complicated than Equation (2.3b). This is because the projectile, which is 
spherical in its own rest frame, appears flattened to an observer at the target, and Equation 
(2.3b) is not a solution of Poisson's equation (second row and second column of Table I) for 
a flattened charge distribution. 

The conditions in the first column of Table I are Lorentz covariant. Thus if ((/? L , A L ) 
satisfy the Lorentz gauge conditions in one Lorentz frame and they are subjected to a 
Lorentz transformation, then the transformed potentials will still satisfy the Lorentz gauge 
conditions. This is not true for potentials satisfying the Coulomb gauge conditions given in 
the second column of Table I. If these potentials are subjected to a Lorentz transformation, 
the transformed potentials will generally not satisfy the Coulomb gauge conditions. In the 
calculation that follows, we begin with a static projectile charge distribution, viewed in the 
projectile rest frame. In this case, the ordinary Coulomb potential satisfies both the Lorentz 
and Coulomb gauge conditions. However when we transform this potential to the target 
rest frame, the resulting potentials satisfy only the Lorentz conditions. If we want projectile 
potentials that satisfy the Coulomb conditions in the target frame, we must find a gauge 
transformation to take us to the Coulomb gauge from the potentials that have been obtained 
by Lorentz transformation from the projectile frame. 

A necessary and sufficient condition that Equations (2.1a,b) and Table I should be com- 
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patible is 

V 2 A(r,t) = ^ L M). (2.6) 

To find a convenient expression for the right-hand side of Eq.(2.6), we start in the rest frame 
of the spherically-symmetric projectile. Let x, y, z(= f ) be position coordinates measured 
relative to the projectile center-of-mass. Then the scalar potential <p(r) satisfies 

V 2 <^(r) = 4vrp p (f) 



whose solution can be expressed [171] as 

1 r r piq-(f-r') 
m = — 2 J dh'J d*q— — Pp (r') 

= - I d 3 q^- r' 2 dr'Uqr')p p (r'). (2.7) 
7T J r Jo 

It is assumed that all the projectile charge is contained within a sphere of radius R p . The 
corresponding vector potential, A(f), is zero. 

If the projectile moves so that its center is located relative to the target by 

r = by + vtz = p + vtz , 

an observer at the target center will measure the scalar potential ip L (r,t) to be 

<p h (x,y,z,t) = jip(x,y,z) (2.8) 

with 

x = x (2.9a) 

y = y-b (2.9b) 

z = j(z-vt) (2.9c) 

7 = (1 - ^r 172 (2-9d) 



Thus we can write 



27 



3 iq-(a;x+(y-b)y+7(2— vt)z) rK 



<p L {x,y,z,t) = — I d 3 q x / r' 2 dr'j Q (qr')p p (r') 



71 J Q 2 Jo 



dqz ^y(z-vt) / d 2 q± x / r ' 2 dr'jo(qr')p p (r') 



'p 



dq z e iqz{z ~ vt) / d 2 



q±- 



giq±-(p-b) 
2 _L f^2 



x j\' 2 dr'jo^q 2 ± + (^r')p p (r'). (2.10) 
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If equation (2.10) is used in equation (2.6), we get 

V 2 A(r,t) = ~J dq z q z e^-^ J A± ^^ 2 r"dr'j + r') p p (r'), 
whose solution can be written 

AM) = ^ [j^ qz{Z ~ Vt) j d2 ^ q ^ ( ^l ?) r' 2 dr'jo (^1 + (^) 2 p p (r') 



2i 



7Tt>C 



giqx-(p-b) 



(gi + (^) 2 )(?i + 



. ^y-y ) ) 



r /2 rfr' Jo (,/gi + (-) 2 r')p p (r') 



7f 



(2.11) 



in which we have replaced the integration variable q z by uo = q z v. 

We can perform the d 2 q± integration in (2.11) by using the relations 

1 1 ~2 



& + (^) 2 <& + 



w\2 



/ 



d^ 



2 iq±-(p-b) 



91 + 



Jo 



7f / \ 7f 



(2.12a) 
(2.12b) 



e iq ± -(p-b) 

£ + (IF J0 



,i + (^) 2 /) = 2^ (M^zM j Jo rn, , . (212c) 



Equations (2.12b,c), which are proven in the Appendix, are valid when r' < |p — b|. This 
condition is satisfied because of our assumption that 6 is large enough so that the projectile 
and target never overlap. Substituting Equations (2.12a,b,c) into Equation (2.11) leads to 



A(r,t) = 4i- 



„ -ILOti — Z 

-e e v 



(2.13) 



x 



Z p e^/M|p-b|\ _ Ko f\u\\p-b\ 



R- 



r f2 dr'j ( « — r' ) p P (r'^ 



An \ 7f y \ t> / j 

This A(r, t) will generate the transformation (2.1a,b) between the Lorentz and the Coulomb 
gauges. 



III. INTERACTION MATRIX ELEMENTS 



The matrix elements of the projectile- charge interaction in the two gauges are given by 



IvHt)] 



fa 



[p T (r)] 7a y? L M) - -[J T (r)] 7a • A L (r,t) 



(3.1a) 
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[V°(t)} 



•yet 



d 3 r 



[p T (r)] 7 ^ c (r,t) - -[J T (r)] 7a .A°(r,t) 



If we use Eqs. (2.1a,b) to express the differences between the potentials, we get 



[vc(t)h*-[vHt)] 



70 



d 3 r 



[P T (r)] 7Q ( -- c ^r) -^[J T (r)] 7 a- VA(r,f) 



(3.1b) 



(3.2) 



Target charge conservation (Equation (1.4)) plus Gauss' theorem, applied to the localized 
target charge density, imply that 



|[J T (r)] 7a -VA(r,t)rf 3 r = - J (V • [J T (r)] 7Q ) A(r, t)d 3 

= iu ia J[p T (r)] ja A(r,t)d 3 r. 



(3.3) 



Thus 



[V c {t)] ia - [VHt)], a = -\j d Z r[p,{v)] ia [^^ + iu^Afof)). (3.4) 

A. Point projectile of charge Zpe. 

If Equations (2.3a, 2.5c, and 1.4) are used in the expression (3.1a) for [V L (t)] 7Q , the result 
can be written in the form 

1 (vt - z) 



[V-(t)\ 



fa 



Z p e J d r J 7a • 



v c 



p-b 



iu ial 2 [(£±)2 + ( vt - z y] 3/2 

Similarly, [V c (t)] 7Q can be obtained by using Equations (2.3b, 2.5a,b and 1.4) in (3.1b): 

1 1 



3/2 



(3.5) 



[V°(t)} 



■ya 



^ I d 3 r J 7a (r) • [ z ( 



+ (p-b)f = 

V (p - b) 2 + (vt - Zf + (vt - 2) y/(p - b) 2 + p - Zf 

(P^y + (vt - z y + {vt - z)^({^y + (vt- z y ) J 



) 



/ 



- Z p e I d 3 r P7 " 

v :p - b) 2 + p - ^) 2 



(3.6) 
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B. Finite spherically-symmetric projectile of charge Zpe. 



If we use A(r, t) from Equation (2.13) in Equation (3.4), we get 



[V°(t)] ia • - [V-(t)U = - f^f^-l) f d 3 r\p T (v)] ia e^ 

^ J — oo \ UJ ' J 



X 



1 Ka 



An 



uj\ \p — b| 



(3.7) 



\u\\p — b| 



r' 2 dr'jo ( i — r' ) p P (r') 



The time-structure of Equation (3.7) suggests that we formulate the expression in terms 
of Fourier transforms 



V c ' l (uj 



> s LP' 



V c ' L (t) = — / due- wt V c ^{u) 
Then Equation (3.7) takes the form 



(3.8a) 



(3.8b) 



[V-{u)\ ia - [V-{u)\ ia = — 



7T (LO. 



TlV V ijj 



l) J d 3 r[p T ( r )], a e 1 ^ 



x 



V^/Mlp-bl^^/MIP-bi 



/' 

Jo 



7"' dr'jo [ i — r' ) p P {r') 



47T \ 7t> / \ v J Jo \ c J 

We see that on-shell (i.e. uj = 0J lct ) interaction matrix elements are the same in the Coulomb 

and Lorentz gauges, which confirms the general result obtained in Section I concerning the 

gauge invariance of onshell interaction matrix elements. 

To find expressions for [V c (u;)] 7a and [V L (c<;)] 7Q , separately, we substitute Equations 

(2. 5c, 2. 10a) into Equation (3.1a), and use Equations (2.12b) and (3.8). A straightforward 

calculation gives 

\P ~ b| 



fa 



2Z e , 
— — / drr 

nv 



p T (r) - 3 [J T (r)] 2 



x e l ^ z K n 



fa 



and then Equation (3.9) gives 



\v^)], a = / dh- 

8n (w- 



[J 



T Jz 



V UJ 



,i — z jy- 



7f 



\p — b| 



(3.10) 



fa 



7f 



x / r' 2 dr'jo fi—r'j p P (r'). 



\uj\ \p — b| 



(3.11) 
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Finally we can proceed, as in Equation (3.6), to express the matrix elements in terms of 
the current density only: 

t^>]- = / • (* + £v) e n> Ko (M^_M) , (3 . 12a) 

2Z p e f 3 r , N1 A ic 2 \ AMIp - b| 



[^ C M] 7 « = / rf 3 r[J T (r)] 7Q • (z + ^v) e'f'ifo ( 



7^ 



vhu ia V u / J \ c 



1 ) / r' (ir'jo ( i — P P ( r ' 



x y rf 3 r[J T (r)] 7Q • Ve*« z ifo (^^— ~J ■ ( 3 - 12b ) 

In the following sections, we will investigate significant differences between expressions 
(3.12a) and (3.12b). 



IV. COMPARISON OF THE STRUCTURES OF V L (w) AND V c {u) 

A. High bombarding energy limits 

At high bombarding energy, where v ~ c, the main bombarding energy dependence enters 
V l (uj) and l /c (u;) via the 7-dependencies exhibited by Equations (3.12a) and (3.12b). Since 

J\^H\ = _ log /MlfLlM) =5 logo, 



7— »oo 

U— >C 



we have 

7 lhn[^H] 7Q = -^^7 /rf 3 r[J T (r)] 7Q -[z + £^V]e^ 
„-> c 2Zpe , / io 



he 2 



log 7 (l - j d 3 r[J T (v)} 7a • ze^ 2 . (4.1) 



Since the <p dependence of [Jt(i")] 7Q • z is given by e t ( M i- M ^)'t > ; the axial symmetry of e l ^ z 
implies that the integral in (|4.1|) vanishes unless M 7 — M a (= //) = 0. Thus we can have a 
log 7 divergence of [V L (co>)] 7a if uj ^ u; 7Q and M 7 = M Q . The effect of this divergence on high 
bombarding energy cross-sections was noted in Reference 
VI below. 



and is illustrated in Section 
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The high-bombarding-energy behavior of [V c (co>)] 7Q , is quite different. The last two lines 
of Equation (3.12b) are independent of 7, and so obviously do not diverge as 7 — > 00. The 
high-7 behavior of the first line is dominated by 

/ d 3 r[Mr)] ja ■ [z + ^V]e^ 

v— >c J 

= -J^ los ~< (i-^)/iMJT«U-^ 

2Z P e log7 /" 3 rT . X1 „ 

The factor of I/7 2 overpowers the log 7 divergence, and the matrix element [V c (u;)] 7Q , is 
seen to be finite at arbitrarily high bombarding energy. 



B. Multipole expansions of V l (uj) and V c (lo) 

It is illuminating to express Equations(3.12a) and (3.12b) in terms of the multipole ex- 
pansion given by WA: 



^ ( M ^ M ) - E ^ (?) E ^ * (?) y f f > < 4 - 2a » 

with Qx^ defined by 



/J 

















X S (27) 2n (n + /i)!n!(A-/i-2n)!' ( ^ 4 ' 2b ' ) 

To expand the third line of Equation (3.12b), we also need the 7 — > 1 limits of Equations 
()4.2a| l4~2b")) . To perform these limits, without affecting the value of v, we allow c — > 00, and 
obtain 



N ' A 



r / — v 2A + 1 v /(A + y u)!(A-/i)! 



. A— u 

u; \ /cur 



A 



x ll7 lv> (4 ' 3) 



If we apply Equations ()4.2a|) and ()4.3|) to Equation (3.12b), we find 

w^)u = ^ E • E [ M<*) + Jtf <"» (?) + x ^ (?) 



M A=|H 



(4.4) 
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where 



IV 



A 2 - /j 2 



Gx + i,» (A + 1) 2 -^ 2 



x 



A y (2A+ 1)(2A- 1) A + l V (2A + l)(2A + 3) 
J rf 3 r[J T ] 7a (r) ■ (V x L)j A (^r)^(f) 



(4.5a) 



X A (M) = x | d 3 r[J T ] 7a (r) ■ Lj A (^r)^(f) (4.5b) 



w 7Q \ 4vr [°° „ 2 ,~ v 
1 I — — / r drpp{r)jo \ir 



47T 



X 



M \ A M 1 /UA A 



7« 



c / V 2A + 1 v /(A + / u)!(A-/x)! 
(^) A x /" rf 3 r[J T (r)] 7a ■ Vr A F M A (f). (4.5c) 



On the other hand, substituting Equation ()4.2aj) into Equation (3.12a) gives 



2Z P e 

hv 



\u\b 

7f 



(4.6) 



ti \=\n\ 

X A (£) and X A (M) in Equation (jOJ) are the same as in Equations ()4.5al and I4.5b|) . but 



X A (G) is defined by 



KiG) = (1-^)^1 d 3 r[J T (r)] 7a ■ Vj 



\ I ~ r } V>). 



(4.7) 



C. Monopole matrix elements 



Inspection of Equation (|4.4j) and (|4.5al I4.5bl and I4.5c|) shows that [V c (u;)] 7a can be 
expressed as 

[V c {u)\ ia = J d 3 r[J T (r)] 7Q • Q c (r) (4.8) 

where Q c (r) is a linear combination of (V x ~L)j\(^r)Y*(r), Lj A (^r)Y^(r), and Vr A Y" / A (f). 
Since 

V-(VxL)j^r)y;(f) = 0, V ■ Lj A (^r)Y A (r) = 0, V • Vr A >; A (f) = 0, 
it follows that 

V-Q c (r) = 0. (4.9) 

Now suppose that <p a and 7 are both states of total angular momentum zero (J = 0). Then 
[Jt(i")] 7Q , is a spherically-symmetric vector field (a central field), and can be written as the 
gradient of a spherically-symmetric scalar field: 



[J T (r)] 7Q = VV(r). 
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Equations (J4~gj) and (JO)) imply that 

[V c (co)} ia = J d 3 rVtp(r) ■ Q c (r) = -J rf 3 r^(r)V ■ Q c (r) = 0. 

so that the matrix elements of V° between any two J = states vanish identically. Thus 
an excited J = state, in a nucleus with a J = ground state, can only be populated 
indirectly, via a multi-step process. However, if [V l (c<j)] 7Q is expressed as 

[V^uo)] ia = j rf 3 r[J T (r)] 7a ■ Q L (r), (4.10) 

then one of the components in the expansion of Q L (r) is Vj\(\u>\r/c)Y^(r) (see Equations 

g2E3D). But 

V ■ V JA (Mr) Y?{t) = V 2 jA (Mr) Yft) = - j x (Mr) 

which is not identically zero. Thus V • Q L (r) is not identically zero and there is no reason 
to expect matrix elements of V L between J = states to vanish. A calculation based on V L 
implies the possibility of one-step population of an excited J = state from a J = ground 
state. 



V. NUMERICAL COMPARISONS OF MATRIX ELEMENTS OF V L AND V c 

The following numerical comparisons refer to Coulomb excitation of the giant dipole 
resonance (GDR) in a 40 Ca target by a 208 Pb projectile. The target transition current density 
matrix elements [Jx(r)] 7Q needed in Equations (4.5) and (4.7) are calculated using Brink's 
model Q| of the GDR, in which unexcited proton and neutron spheres undergo harmonic 
oscillations relative to each other. The radial densities of the spheres are obtained by filling 
the lowest available shell-model single-particle states. The relative harmonic oscillations of 
the proton and neutrons spheres are characterized by the eigenstates (p^f(r pn ) of a three- 
dimensional harmonic oscillator. The ground state is . The degenerate first excited 
states that can be populated by Coulomb excitation are ^j} 1 and + <p°l\) /y/2, the 
combination of 0^ symmetric under reflection across the reaction plane. Further details on 
he calculation of matrix elements using the Brink model eigenstates are given in reference 

In the calculations represented in Figure 1, it has been assumed that the 82 protons of 
208 Pb are distributed uniformly within a sphere of radius 7.5 fm (as seen by an observer 



16 



moving with the projectile). The parameters used for the calculations imply that b = 20 fm 
and u ia = 11.7 MeV/h. It can be verified by inspection of Figure 1 that, in every case, V c 
and V L agree at uj = 11.7 MeV/K, as required by Equation (3.9). 

In Figure 1, /x = refers to the transition 0q° — > 0q 1 and fi = 1 refers to the transition 
0o'° — > + 4>li) /y/2- In every case, the solid line gives V c (uj), whereas the dashed line 
gives V L (u). 

The first observation concerning Figure 1 is that at a bombarding energy of E/A=100 
MeV, there is little difference between V c and V L . However, at E/A=10 GeV, the difference 
is pronounced. The most striking difference is the very strong increase with bombarding 
energy of V L for the /i = transition. This is an expression of the logarithmic divergence 
(with increasing 7) of the /x = matrix elements of V L , as shown in Equation f)4.1j) . It is 
clear that, at high bombarding energy, calculations using V L will ascribe a much higher role 
to fi = transitions than will calculations using V c . 

Figure 1 also shows that, at high bombarding energy, the use of the Coulomb gauge 
leads to an interaction that is more adiabatic than predicted by the Lorentz gauge. The 
wider spread of V L (u), as a function of u>, compared to the spread of V c (u), shows that 
the impulse in the Lorentz gauge is sharper than the Coulomb gauge impulse. For example, 
Figure 2 shows the same comparison as in the E/A = 10 GeV, /i = 1 plot of Figure 1, but 
in the time domain. The sharpness of the pulse provided by V L (t), compared with that 
provided by V c (t), is evident. Note that in this particular case, [V L (— uj)] ya = [V L (u)] ia 
which has the consequence that [V L (t)] 7Q is real. However, this is not true of [V c (c<j)] 7Q , and 
thus [V c (£)] 7 a has both real and imaginary parts. 

VI. GAUGE-DEPENDENT EFFECTS IN COUPLED-CHANNEL 

CALCULATIONS. 

The most significant test of the differences between V L and V c is in the calculation of RCE 
cross-sections, since the cross-section is the point where theory and experiment intersect. In 
this Section, we compare cross-sections calculated with these two interactions, for the 208 Pb 
+ 40 Ca system described in Section V. We have performed coupled-channel integrations of 
the time-dependent Schrodinger equation. The target states included in the calculation span 
the 0-phonon, 1-phonon, and 2-phonon states of the 40 Ca GDR. The methods used to do 
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the numerical Fourier transform needed in Equation (3.8b), and to integrate the coupled 
equations, are described in Reference Q. We also describe there the integrations over 
impact parameter needed to calculate the cross-section. 

Figure 3 compares the calculated cross-sections for the population of the six reflection- 
symmetric 1- and 2-phonon states that can be reached via RCE, as functions of the kinetic 
energy of the 208 Pb projectile nucleus. In every case, a solid line is used to show the result of 
the calculation using V c (t), and a dashed line is used to show the result of the corresponding 
calculation using V L (t). 

For the one-phonon states, the situation is similar to that shown in Figure 1. The 
two sets of calculations agree at low bombarding energy (E/A < 1 GeV). However, at 
high bombarding energy (E/A > 5 GeV), V L (t) predicts much greater cross-sections for the 
population of the J = 1, M — state than does V c (t). The situation is more complicated for 
the two-phonon states. The cross-sections are much smaller than for the one-phonon states, 
and depend upon multiple excitation processes. Also, the truncation of our calculation at 
two phonons introduces an additional element of uncertainty into our two-phonon cross- 
sections. However, it is noteworthy that the 2-phonon \M\ = 1 state is also strongly favored 
by V L (t) at high bombarding energy, compared to V c (t). The reason is that this state is 
mostly populated in two-step processes, such as 

(J = M = 0) -»• (J = 1, M = 0) -»• (J = 2, |M| = 1) 

and 

(J = M = 0) -»• (J = 1, |M| = 1) -»• (J = 2, \M\ = 1). 

In both AM = transition is involved, and it is the logarithmic divergence (with 

increasing 7) of the AM = transition amplitude that leads to the strongly increasing 
cross-section when V L (t) is used. But the J = 2, \M\ = 2 state is reached mainly by 

(J = M = 0) -> (J = 1, \M\ = 1) -> (J = 2, \M\ = 2), 

in which there is no AM = step, and so the J = 2, \M\ = 2 state is not strongly favored 
by V L (t) at high bombarding energy. Thus we see that for E/A > 5 GeV, V c {t) and V L (t) 
predict very different cross- sect ions when used in truncated coupled-channel analyses. A 
calculation using V L (t) will predict stongly enhanced cross-sections for populating any state 
that can be reached via a one-step or two-step process involving a AM = transition. 
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The stongest transition illustrated in Figure 3 populates the one-phonon \M\ = 1 state. It 
is predominantly a one-step transition, and so is well described by the Born approximation. 
According to Equation f)4.4j) or (|4.fi|) . this implies that the 6-dependence of the transition 
probability is given by 



whose the integral over impact parameter diverges in the 7 — > 00 limit (cf Equations (3.1), 
(3.6) of Reference Q].) The physical reason for this divergence is that as 7 increases and the 
electromagnetic pulse becomes more strongly retarded, it becomes flatter and its influence 
extends for a longer distance away from the trajectory of the projectile. Thus, in the integral 
over b, larger values of b play a more important role as 7 increases, and in the 7 — > 00 limit, 
the 6-integral diverges. This occurs only for AM = ±1, because the flat pulse is spatially 
axially symmetric (AM = 0) and the intrinsic spin of the photon transfers AM = ±1. The 
one-phonon \M\ = 1 cross-section shown in Figure 3 exhibits this logarithmic-like increase 
with bombarding energy. 

However, the log 7 divergence of the off-shell /i = AM = matrix element, discussed 
in Section IV. C, is a much more serious divergence. It occurs for every b, and thus the 
6-integrated cross-section would be a divergence of still higher order. This is the diver- 
gence that is introduced when the Lorentz gauge is used in a coupled-channel calculation, a 
divergence which is not present when the Coulomb gauge is used. 

The expression for V L (u) given by Equation (3.12a) depends upon the total projectile 
charge, but not on the radial dependence of this charge. However there is a dependence on 
the projectile radial charge density in V c (oj) given by Equation (3.12b). Fortunately, this 
dependence is very weak. The calculations using V c (uj) whose results are shown in Figure 
3, were done using a radial charge distribution which was constant from the center out to 
7.5 fm. By way of comparison, we have repeated these calculations, assuming that all the 
charge of the 208 Pb nucleus is concentrated at its center. In all cases, we found that the 
calculated cross-sections changed by no more than a few tenths of a percent. As might be 
expected, the transition amplitudes are more sensitive to projectile radial density at small 
impact parameters. However, when the integration over all impact parameters is done, the 
residual effect on the cross-section of changes in the projectile radial charge density is very 
small. 
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VII. CONCLUSIONS AND DISCUSSION 

We have studied the relationship between the electromagnetic potentials in Lorentz gauge 
and Coulomb gauge for the fields encountered in relativistic Coulomb excitation, and have 
found expressions in the two gauges for the interaction between projectile and target. At 
bombarding energies above about 2 GeV per nucleon, we have found significant differences in 
excitation cross-sections when the two gauges are used in coupled-channel time-dependent 
calculations, especially for processes involving AM = transitions. Since there is no a 
priori reason for using one gauge rather than another, this lack of gauge invariance reveals 
a weakness of the coupled-channel time- dependent approach to relativistic Coulomb exci- 
tation. It demonstrates that the common practice of relying totally on calculations in the 
Lorentz gauge cannot be justified. 

We have discussed one cause of lack of gauge invariance, the truncation of the set of target 
states used in the coupled-channel calculation. Another cause is revealed by consideration of 
the Hamiltonian for a target charged particle moving in the fields (v?(r, t), A(r, t)) produced 
by the projectile: 

H = J- f p _ ^fcl^y + & p(r,t) + C/ nuc (r) 



2m \ c 
= ^p 2 -^(p-A(r,t) + A(r,t).p) + ev ,(r,t) 

+ _^A(r,t)-A(r,t) + C/ nuc (r) (7.1) 



The terms in Equation (7.1) that are linear in e, 



e 



p • A(r, t) + A(r, t)-p ) + ey9(r, t), 



2mc 

give rise to the interaction (1.3) that has been the basis of this study. However the term 
quadratic in e, 

^A(r, t ).A(r, t ), 

has not been included in our calculation, or in other studies published to date on relativistic 
Coulomb excitation. But the equation 

e(A(r,t) + VA(r,t))\ 2 i£ a(m) f eA(r,t)\ 2 

' e ^ ifj(r, t) = e ^ p — - ip(r,t), 



c J \ c 

which is an essential component in the argument about the gauge invariance of the observable 
consequences of the Schrodinger equation, relies for its validity on the presence of the A ■ A 
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term. Thus it is clear that no gauge invariant theory of relativistic Coulomb excitation can 
be constructed without inclusion of this term. We will address this issue in a forthcoming 
publication. 



APPENDIX: TWO INTEGRAL RELATIONS 



Define 1(f) by 



1(f) EE / £^ e^<P^ jo ( f^ql + A| ) . (A.l) 



Operate on this expression with V 



2. 



+ A 



(q ± -p+A 2 i) 



- {A - A?) / ^| ^^- b) Jo [ry/& + A| ) . (A.2) 
In the first term of Equation (|A.2|) . replace j ^f -y/g^ + A|J by 

Jo (ryj<j[+ A|J = ^ J sin ^0#e 
where p ee f sin#(cos0x + sin0y) and 5 = fcosO. Then Equation (jA.2|) becomes 

V?J(f) = --^ / sm6d9d<j) J d\ ± e^iP-^'P) e ~ iX2 ~ z 

VpJ(f) + (Ag - A?)J(f) = --^2 / sin 6d6d0e- tX2i 5 (p - b - p). (A.3) 

47r y 

In our application, |p| is bounded by the radius of the target and |p| is bounded by the 
radius of the projectile. The condition that the projectile and target do not overlap implies 
that |p — b| > |p| for every orientation 6, (ft, and so p — b — p is never zero. Thus the 
^-function on the right-hand side of Equation ()A.3j) vanishes, and we get 

V?J(f) + (A* - A?)I(r) = 0. (A.4) 
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We have two applications of Equation (|A.4|) . In one of them, Af = A| = {uj/{jv)) 2 . In 
the other, = (u/v) 2 , X 2 = {uj/{^v)) 2 . Thus in both cases, (X 2 , — Xf) < 0, and we can 
write the general solution of Equation (|A.4|) as 




(A.5) 



where a and (3 are independent of r. The first term in Equation (|A.5|) is finite at r = 0, 
whereas the second term diverges there. But if we set f = in the definition (|A.1J) . we get 

J (°) = / ^~TT2 e^<P-V = 2nK (\X 1 \\p-b\) (A.6) 
which is finite in our \p — b| > situation. Thus (3 in Equation ()A.5|) vanishes and 

I(0)=aj (0)=a = 27iK (\X 1 \\p-b\), 

leading to 

I(r) = j e^P'V Jo (r^ql + xij=27rK (\X 1 \\p-b\)j (iryj\* - A|) . 

9± 1 (A.7) 
Thus we get the two special cases: 

/ ^w- ^ * - (^) < A -> 

/ ^ ^> * W*^) ^ (^) * B 1 ) ■ (A*) 
These equations are valid when f < |p — b|. 
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Figure Captions 

Figure 1. Plots of matrix elements of V c (u) (solid lines) and V l (uj) (dashed lines) 
connecting the Ca ground state to the one-phonon giant dipole resonance states. The 
projectile is a 208 Pb nucleus with bombarding energies per nucleon specified within the 
figure frames, n — corresponds to the M = one-phonon state; \x = 1 to the reflection- 
symmetric \M\ — 1 state. 

Figure 2. Comparison of V c (t) (solid line) and V L (t) (dashed line), corresponding to the 
E/A = 10 GeV, fi = 1 example of Figure 1. 

Figure 3. Calculated coupled-channel Coulomb excitation cross-sections for one- and 
two-phonon states of the GDR in 40 Ca with 208 Pb projectiles. The solid curves correspond 
to calculations in which the Coulomb gauge has been used, the dashed curves to calculations 
with the Lorentz gauge. 
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Figure 1 Bayman and Zardi Gauge dependence. 
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